Classical simulations including electron correlations for sequential double ionization 
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With a classical ensemble model that including electron correlations during the whole ionization 
process, we investigated strong-field sequential double ionization of Ar by elliptically polarized pulses 
at the quantitative level. The experimentally observed intensity-dependent three-band or four-band 
structures in the ion momentum distributions are well reproduced with this classical model. More 
importantly, the experimentally measured ionization time of the second electrons [A. N. Pfeiffer et 
al., Nature Phys. 7, 428 (2011)], which can not be predicted by the standard independent-electron 
model, is quantitatively reproduced by this fully classical correlated model. The success of our work 
encourages classical description and interpretation of the complex multi-electron effects in strong 
field ionization where nonperturbative quantum approaches are currently not feasible. 

PACS numbers: 32.80.Rm, 31.90.+S, 32.80.Fb 
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Among various intense laser-induced phenomena, 
strong field double ionization (DI) is one of the most 
important and fundamental processes. During the past 
decades, a great number of experimental as well as theo- 
retical studies have been performed on this area. It has 
been known that DI proceeds either sequentially or non- 
sequentially. In nonsequential double ionization (NSDI), 
the second electron is ionized by the recollision of the 
first tunneled electron [l| . Because of this recollision, the 
two electrons from NSDI exhibit a highly correlated be- 
havior 043- I n sequential double ionization (SDI), it is 
usually assumed that no correlation exists between the 
two electrons, and thus the ionization of the electrons 
can be treated as two independent tunneling-ionization 
steps. However, this assumption has been called into 
doubt by recent experiments @, [1] ■ In Ref . 7] , it has 
been shown that there is a clear angular correlation be- 
tween the two electrons from SDI, which implies that the 
successive ionization steps are not independent in SDI. 
In Ref. [H, the authors found that the ionization time 
of the second electron from SDI is much earlier than the 
prediction of the independent-electron model. These ob- 
servations declare that the electron correlations in SDI 
should be reexamined carefully. 

Theoretically, an accurate description of the electron 
correlations in DI needs full quantum theory. However, 
due to the enormous computational demand, the solu- 
tion of the time-dependent schrodinger equation of the 
two-electron system in the strong field, especially in the 
case of elliptically polarized laser fields, is not feasible 
currently. Instead, a fully classical treatment of the two- 
and multi-electron systems proposed by Eberly et al. has 
been well established [§-[l2[ . During the past decade, this 
model has been successful in exploring the strong-field 
ionization processes at the qualitative level. However, 
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it fails when made a quantitative comparison with ex- 
periments. For example, in ref. fl3| it has been shown 
that the saturation intensity for SDI of Ar in the clas- 
sical simulation is much higher than the experimental 
data [1] , and the amplitude of oscillation in the ratio of 
the parallel to antiparallel emitted SDI counts as a func- 
tion of inte nsity is much larger than the experimental 
observation [1J|. It remains a question: can a classical 
treatment describe the strong field processes at the quan- 
titative level [TH? In this Letter, we gave an affirmative 
answer. With a delicate modification to Eberly's model, 
we performed a first quantitative simulation on strong 
field SDI of Ar in the close-to-circular laser fields. Our 
numerical results exhibit intensity-dependent three-band 
or four-band structures in the ion momentum distribu- 
tions, which are consistent with the recent experimental 
results @. Especially, the experimentally measured ion- 
ization time of the second electron @ , which can not be 
predicted by the standard independent-electron model, 
is reproduced surprisingly well with this purely classical 
model which fully takes into account the electron corre- 
lations during the entire SDI process. The quantitative 
agreement between our numerical results and experimen- 
tal data indicates that a classical treatment is a good ap- 
proximation in describing the strong-field SDI and has 
the potential to shed light on the subtle multi-electron 
effect in strong-field double and multiple ionizations. 

In the classical picture of a two-electron system, one 
electron often drops deeply into the nuclear potential 
well, leading to the autoionization of the other electron. 
In Eberly's classical model [t| [l(| EH j a soft-core poten- 
tial is introduced for the ion-electron interaction to avoid 
autoionization. However, in this soft-core potential clas- 
sical model (SPCM), the first and second ionization po- 
tentials can not be matched with those of the investigated 
target. In fact, the first electron often ionizes leaving the 
second electron with an energy much lower than the sec- 
ond ionization potential of the target [l7|, EH ■ In tunnel- 
ing and over-the-barrier ionization, the ionization rate is 
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very sensitive to the ionization energy. Thus, the SPCM 
may be deficient in describing strong field double ion- 
ization that occurs through tunneling or over-the-barrier 
escape. However, this deficiency can be avoided with the 
Heisenberg-core potential, which not only prevents au- 
toionization but also gives the ground-configuration ener- 
gies of the multi-electron atoms [19J. In the past decades, 
this potential has been extensively employed in the clas- 
sical investigations of atomic and molecular collisions and 
laser- matter interactions [13, HH ■ Here, we employ the 
Heisenberg-core potential instead of the soft-core poten- 
tial in the classical two-electron model to study SDI of 
Ar by the elliptical laser pulses. 

The Hamiltonian of the two-electron atom in the 
Heisenberg-core potential classical model (HPCM) is 
(atomic units are used throughout this Letter until stated 
otherwise): 

"o = ^ + £h^ + f + M^)] (i) 

where and Pi are the position and canonical momen- 
tum of the ith electron, respectively. Vh Wi ,Pi ) is the 
Heisenberg-core potential, which is expressed as [191 ]: 

V H (n, Pt ) = j^exp{a[l - (^p) 4 ]} (2) 

The parameter a indicates the rigidity of the Heisenberg 
core and is chosen to be 2. For a given a, the parameter 
£ is chosen to match the second ionization potential of 
the target, i.e., it is set to make the minimum of the one- 
electron Hamiltonian [Hi = 77 + ^ + Vh(ti,Pi)} equal 
to the second ionization potential of Ar (-1.01 a.u.). For 
a— 2 we obtain £=1.225. For the two-electron atom, in 
refs. jl9l - [2l1 |. the ground-state energy of the system is 
obtained by minimizing the Hamiltonian Hq. In that 
configuration, the two electrons locate at opposite sides 
of the nucleus, and they are stationary to each other 
[lijj . Differently, in our calculation, we employed the 
approach used by Eberly et al. [16| to determine the 
ground-state energy of the atom, i.e., by inputting an 
energy that equals to the sum of the first and the second 
ionization potentials of the target (-1.59 a.u. for Ar). 
The initial distributions of the ground-state atom in the 
phase space are obtained with the approach in SPCM 
US EH- Figures. 1(a) and 1(b) display the position dis- 
tributions of the ground-state atom. We state that in 
our calculation, the second ionization potential is related 
to the minimum of Hi while the ground-state energy of 
the atom is determined by the "input" energy, which is 
higher than the minimal value of Hq. Thus, the two elec- 
trons are not fixed at two certain points but distributed 
in a finite region of the phase space [see figs. 1(a) and 
1(b)]. This relaxation of the position of the electron pairs 
in the phase space allows sufficient electron correlations 
in the initial state and during the ionization of the first 
electron. 




Energy (a.u.) Time (Is) 

FIG. 1: (color online) (a)(b) The initial position distributions 
of electrons in the HPCM, where the initial energy of the sys- 
tem is set to be -1.59 a.u. The parameters of the Heisenberg 
core are a — 2 and £ = 1.225. Here (xi, yt, Zi) represent the 
coordinates of the ith electron in the directions of the x-, y-, 
and z-axis, respectively, (c) The ionization energy distribu- 
tions of the first (the green curve) and the second (the red 
curve) electrons, (d) The energy evolution of the electrons 
for a typical SDI trajectory. 



The Hamiltonian of the two-electrons in the presence 
of the laser field is 

H = H + (r 1 + r 2 ) . E(t) (3) 

where E(t) is the electric field of the laser pulses. The 
evolution of the system in the laser field is determined by 
the following equations: 

dri dH dpi dH 
dt dpi ' dt dri 

The electric field is given as E(f) = f(t)[ .^ cos(u)t + 

l) x + ^ +1 s » n M + l)yh where i{t)=E exp{-±(±:) 2 } is 
the field envelope, w, e and <p are the laser frequency, 
the ellipticity and carrier-envelope phase (CEP), respec- 
tively. 2v7n2r denotes the pulse duration (FWHM). 

First of all, we perform a calculation to test how well 
the first and second ionization potentials of Ar have been 
produced by this HPCM. The two-electron atoms are ex- 
posed to a 33-fs, 788-nm laser pulse with the intensity 
1=4.0 PW/cm 2 and ellipticity e = 0.77. We trace the en- 
ergy evolutions of the two electrons and record the energy 
of the second electron when the first electron is ionized. 
This energy is assumed to be the second ionization po- 
tential of the model atom, which is shown in fig. 1(c). 
The first ionization potential of the model atom is also 
shown in fig. 1(c), obtained by subtracting the second 
ionization potential from the initial energy of the two- 
electron atom (-1.59 a.u.). It is clearly shown that the 
first and second ionization potentials of the model atom 
are about -0.57 a.u. and -1.02 a.u. (almost at the bot- 
tom of Hamiltonian -Hi), respectively, very close to the 
realistic first and second ionization potentials of Ar. 
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FIG. 2: (color online) The ratio of SDI electron counts of 
parallel and antiparallel emissions along the minor elliptical 
axis (x axis) of the laser pulses, as a function of laser intensity. 
The laser parameters are the same as those in ref. [IJ], i.e., 
the duration, wavelength and ellipticity of the pulses are 7 fs, 
740 nm and 0.78, respectively. The results are obtained by 
averaging over the CEP between and 2ir. 



A recent experiment has shown that in SDI of Ar by 
strong elliptical laser pulses the ratio of the parallel and 
antiparallel electron emissions along the minor elliptical 
axis exhibits an oscillating behavior as a function of laser 
intensity |14j . This behavior has been predicted by the 
SPCM and explained as a multi-electron effect beyond 
independent-electron assumption [l3j. We display our 
calculations with the HPCM in fig. 2, where the oscil- 
lating behavior is also clear. It implies that the electron 
correlations in SDI are included and represented in this 
HPCM. One can clearly see from fig. 2 that at the high 
laser intensities the oscillating curve is a bit below 1, 
meaning the two electrons prefer to emit into the op- 
posite hemispheres. This behavior wel l ag rees with the 
experimental data (see fig. 2(a) of ref. fl4|). 

In figs. 3(a) and 3(d) we display the ion momen- 
tum distributions in the polarization plane for SDI of Ar 
by 33-fs laser pulses with laser intensities 1.0 PW/cm 2 
and 4.0 PW/cm 2 , respectively. The ellipticity and wave- 
length are respectively 0.77 and 788 nm, same as those in 
the recent experiment Q. In the direction of the minor 
elliptical axis {x axis), the distribution exhibits a three- 
band structure at the relatively low laser intensity while a 
four-band structure at the relatively high laser intensity. 
In the direction of the major elliptical axis {y axis), as 
shown in figs. 3(c) and 3(f), the spectra show a gaussian 
shape for both intensities. These results agree well with 
the experiment Q . The origin of the four-band structure 
is interpreted as due to the different values of the electric 
field at which the two electrons are released 0, [23[ ■ Be- 
cause of the ellipticity, the electrons preferentially ionize 
along the major axis, leading to the electrons with final 
momenta along the minor axis. The two outer bands in 
fig. 3(d) correspond to the events where the two electrons 
emit into the parallel directions whereas the two inner 
bands result from the events where the two electrons re- 
lease into the antiparallel directions. At the relatively 
low laser intensity, the momentum amplitudes of the two 




FIG. 3: (color online) The ion momentum distributions in 
the polarization plane for SDI of Ar at laser intensities of (a) 
1.0 PW/cm 2 and (d) 4.0 PW/cm 2 . (b) and (c) show the ion 
momentum spectra along the minor and major axis of the 
polarization, obtained by integrating the distributions in (a) 
over p y and p x , respectively, (e) and (f) Same as (b) and (c) 
but for the distribution in (d). 



electrons are almost the same because both electrons are 
ionized around the pulse centre [1] . Thus the antiparallel 
electron emissions result in the nearly zero momentum 
of ion, leading to the three-band structure in ion mo- 
mentum distribution [fig. 3(a)]. The above picture for 
the intensity-dependent ion momentum spectra has been 
confirmed by back analyzing the SDI trajectories of our 
calculations (not shown here). 

Figure. 4(a) displays the correlated radial momen- 
tum distribution of the electron pairs from SDI of Ar. 
The laser parameters are the same as those in fig. 3(d). 
Note that in this calculation the focal volume effect as- 
suming a Gaussian beam profile has been considered. 
The momentum p' r in fig. 4(a) is defined as p' r — 



l)Py, which is an injective func- 



tion of time under the condition that the electron ionized 
before the peak of the pulse @. The repulsion behavior 
along the diagonal means that the two electrons achieve 
different final momenta at the end of the pulses, implying 
the different release times of the two electrons. In ref. 0], 
the ionization times of the two electrons in SDI are read 
from the electron's final momentum p' r . It is found that 
the ionization time for the first electron agrees well with 
the prediction of the independent-electron model. How- 
ever, the ionization of the second electrons occurs much 
earlier than the prediction of the independent-electron 
model. With the same procedure as that in ref. [8], we 
extract the release times of two electrons in our classical 
simulations where the electron correlations are included 
during the entire process of SDI. The results are shown 
in fig. 4(b), where the experimental data from ref. Q are 
also displayed. Note that for our numerical results, the 
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FIG. 4: (color online) (a) The electron correlation spectrum 
for the momentum p' r . (b) The release times of the first (blue 
squares) and the second (red circles) electrons in SDI, which 
are extracted from the electron momentum p' r (see text for 
detail). The experimental data from ref. 0] (triangles) are 
also shown in (b). For comparison, we also displayed the 
numerical results from the SPCM (the magenta circles and 
the cyan squares). 



laser intensities in fig. 4(b) are scaled with a constant fac- 
tor of 0.82, which is well in the experimental uncertainty 
range. Surprisingly, the release times from our classical 
simulations agree excellently well with the experimental 
data, both for the first and the second electrons. For com- 
parison, we have repeated the calculations above with 
the SPCM. The ionization times of the first and the sec- 
ond electrons from SPCM [as shown in fig. 4(b)] devi- 
ate seriously from the experimental data. Here, we only 
displayed the results at intensities above 4.0 PW/cm 2 
because in the SPCM the DI yield is very low at rela- 
tively low intensities (the DI probability does not exceed 
1% until the laser intensity reaches 4.0 PW/cm 2 ). Dur- 
ing the past decades, classical methods have been widely 
employed to describe the strong field phenomena. How- 
ever, they are standing at the qualitative point of view. 
As far as we know, our calculation is the first successful 
example of quantitative investigation of the strong field 
DI with a classical treatment. 

For the theoretical treatment of the independent- 
electron model in ref. Q, the Coulomb interaction be- 
tween the ion and the escaping electron is neglected. Here 
we estimate how strong the Coulomb interaction affects 
the release-time reading. It has been demonstrated that 
the angular distributions of the electrons can serve as 
a signal to estimate the importance of the Coulomb in- 
teraction on the electrons @, [HJ]. In fig. 5 we display 
the angular spectra of the first and the second electrons. 
Note that the angular spectra peak at 180° (and 0°) if 
there is no interaction between the ion and the escap- 
ing electron. Fig. 5 shows that the spectrum of the 
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FIG. 5: (color online) The angular distributions of the two 
electrons from SDI by the 33- fs elliptical (e=0.77) laser pulses. 
The solid lines correspond to the Gaussian fitting of the dis- 
tributions. 



first electron (the red curves) peaks at an angle slightly 
deviating from 180°, indicating a weak Coulomb correc- 
tion on the electron trajectory [24| [25|. For the second 
electron (the green curves), however, the spectrum peaks 
almost at 180°, implying negligible Coulomb interaction 
between the ion and escaping electron. These results in- 
dicate that the influence of the Coulomb attraction on 
the second electron is even weaker than that on the first 
electron. Thus, it confirms that the Coulomb attrac- 
tion has negligible contribution to the deviation between 
the experimental measurement and independent-electron 
prediction for the ionization time of the second electron 
reported in Ref. 

It has been proposed that inelastic tunneling [26|, [53] , 
where the first electron escapes leaving the core in an 
excited state, is important in double and multiple ioniza- 
tion at high laser intensities. However, the most typical 
picture of the ionization process revealed by our classical 
model is that the first electron ionizes leaving the second 
electron at the ground sate of Ar + , as shown in fig. 1(d) 
[see the inset of fig. 1(d)]. Note that the second electron 
stays at a state with an energy higher than -1.01 a.u. af- 
ter the first ionization for a small part of the SDI trajecto- 
ries shown in figure 1(c). This "excitation" is responsible 
for the oscillation behavior of the parallel-antiparallel ra- 
tio (see fig. 2) [l3| . However, it is not responsible for the 
early release time of the second electron because the sec- 
ond ionization time shown in fig. 4(b) does not change 
when these trajectories are excluded. Another reason for 
the deviation of the independent-electron prediction of 
the second ionization from the measured data [8| possi- 
bly ascribes to the empirical formula for ionization rate 
[281 ] . which may be inaccurate in the experimental con- 
dition where nonadiabatical effect is important to the 
ionization process [2{|. Then, it is not clear to what ex- 
tent the deviation of the second ionization time comes 
from the inaccuracy of the ionization rate formula and to 
what extent the electron correlations influence the ion- 
ization of the second electron. These questions call for 
more detailed investigations. 

In conclusion, we have investigated SDI of Ar by the el- 
liptical laser pulses with the HPCM in which the electron 
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correlation is included during the entire process. The 
experimental observed ion momentum spectra and oscil- 
lating behavior of the ratio of the antiparallel to paral- 
lel electron emissions are well reproduced by this classi- 
cal model. Especially, the measured ionization time of 
the second electron, which strongly deviates from predic- 
tion of the standard independent-electron model, is excel- 
lently reproduced by the HPCM. The quantitative agree- 
ment between our classical calculations and experimental 
results provides strong supports to the classical treat- 



ment of the multi-electron processes induced by strong 
laser fields, which is currently indispensable because the 
nonperturbative quantum treatments of the complex ef- 
fect are not feasible. 
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